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Abstract 

Resonant x-ray reflectivity measurements from the surface of liquid Bi22ln78 find only a modest 
surface Bi enhancement, with 35 at% Bi in the first atomic layer. This is in contrast to the 
Gibbs adsorption in all liquid alloys studied to date, which show surface segregation of a complete 
monolayer of the low surface tension component. This suggests that surface adsorption in Bi-In is 
dominated by attractive interactions that increase the number of Bi-In neighbors at the surface. 
These are the first measurements in which resonant x-ray scattering has been used to quantify 
compositional changes induced at a liquid alloy surface. 

PACS numbers: 61.25.Mv, 68.10.-m, 61.10.-i 



1 



I. INTRODUCTION 



Current treatments of the thermodynamics of surface phenomena in solutions rely heavily 
on the original works by Gibbs in 1878, and one of the most familiar corollaries is the Gibbs 
adsorption rule. In its simplest invocation, the Gibbs rule states that in a binary liquid, 
the species having the lower surface tension will segregate preferentially at the surface. 
This apparent simplicity is deceptive: a survey of the literature reveals a hundred years' 
debate over the application of the Gibbs adsorption rule,~ not to mention its extension to 
multicomponent systems^ and crystalline surfaces,— and its connection to atomistic models.— 

Experimental investigations of the validity of the Gibbs rule encompass measurements of 
adsorption isotherms,— surface tension^^ and surface composition^ in a variety of systems. 
Unfortunately, many of the liquids studied are too complicated for the simplest formulations 
of the Gibbs rule. Liquid metal alloys are in many ways ideal for such studies. Miscible alloys 
exist which behave as ideal liquids, while in other systems strongly attractive or repulsive 
heteroatomic interactions can be studied. Perhaps an even more important advantage of 
liquid metals is that the compositionally inhomogeneous region at the surface is known in 
some cases to be confined to an atomic layer. This is commonly assumed in calculations of 
Gibbs adsorption that take a model of a physical surface as their starting point. 

For example, x-ray reflectivity, ion scattering, and Auger electron spectroscopy measure- 
ments of liquid Ga84lni6 found a 94% In surface monolayer, as expected given this alloy's 
positive heat of mixing.^^ Subsequent layers have the bulk composition. Similar studies of di- 
lute liquid Bi-Ga (< 0.2 at% Bi) likewise found surface segregation of a pure Bi monolayer.— 
Even when the repulsive interactions between Ga and Bi cause more Bi-rich alloys to undergo 
additional phase separation above 220°C, where a 65 A thick inhomogeneous Bi-rich region 
forms, the pure Bi surface monolayer persists . -^^i-^^ In Ga-Bi, then, repulsive heteroatomic in- 
teractions substantially change the surface composition profile, but do not defeat the Gibbs 
adsorption. 

The effect of attractive heteroatomic interactions remains an open question. In alloys 
such as Bi-In, attractive forces between the two species produce a number of compositionally 
ordered phases in the bulk solid. It is therefore conceivable that Bi-In pairing may exist at 
the liquid surface, and compete with surface segregation. Our recent temperature dependent 
x-ray reflectivity measurements of liquid Bi-In alloys having 22, 33, and 50 at% Bi revealed 
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structural features not found in elemental metals or in the Ga alloys discussed above.— 
As we will show, those data were suggestive of Bi-In pair formation along the surface- 
normal direction. However, since the technique did not measure the Bi surface concentration 
directly, other interpretations of the data were also possible. 

II. EXPERIMENTAL SETUP 

A complete characterization of surface composition requires both elemental specificity and 
A-scale structural resolution along the surface-normal direction, which is difficult to achieve 
experimentally. Auger electron spectroscopy, which satisfies the first of these requirements, 
is hampered by contributions from the bulk liquid.— X-ray refiectivity by contrast is a surface- 
sensitive probe. In the kinematic limi1?i^ the refiected intensity, measured as a function of 
momentum transfer normal to the surface, is proportional to the Fresnel refiectivity Rp 
of a homogeneous surface-^ 



Here pefr represents an effective electron scattering amplitude that combines the electron 
density profile with the scattering form factor, and poo is the density of the bulk. The electron 
density variations that produce modulations in the refiectivity may result from changes in 
either the composition or the mass density. Thus, inference of surface composition from the 
measured refiectivity is sometimes ambiguous. 

III. RESONANT X-RAY SCATTERING 

This disadvantage can be overcome with the application of resonant x-ray scattering. The 
effective electron density of a scattering atom depends on the scattering form factor f{q) + 
f'{q, E) ^ Z + f\E). When the x-ray energy is tuned to an absorption edge of a scattering 
atom, the magnitude of /' becomes appreciable, producing changes in contrast between 
unlike atoms.— With one exception,— resonant x-ray scattering measurements reported in 
the past have been confined to studies of solids and bulk liquids, due to the difficulty of the 
experiment. The present report is the first to find compositional changes induced at a liquid 
surface. 
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FIG. 1: (a) Energy scan in transmission through Bi foil, (b) Fresnel-normaUzed x-ray reflectivity 

of Bi22ln78. o: 12.5 keV (two independent measurements); •: 13.4 keV; ( — ): 12.5 keV fit; ( ) 

13.4 keV fit. (c) Best fit real space number density profile relative to bulk atomic percent. ( — ): 

Total number density; ( ): In density; ( ): Bi density, (d) Surface energy as a function of 

surface Bi concentration Xo.Bi for the bulk concentration XBi = 0.22, according to Eq. ^ ( ): 

w = . {—): w = -WksT. 

IV. SAMPLE PREPARATION 



The molten Bi22ln78 sample was maintained at T = 101°C and P = 5 x 10"^*' Torr 
within an ultra high vacuum chamber, and periodically sputter cleaned with Ar"*" ions. 
Reflectivity measurements were performed at beamline X25 at the National Synchrotron 
Light Source. The spectrometer has been described previously,"^^. except that here a double 
Si(lll) crystal monochromator was used to provide an energy resolution of 9 eV. For the 
present work, we compare reflectivity measured at 12.5 keV to measurements made at the Bi 
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Ljjj edge at 13.421 keV. The energy was calibrated by transmission through a Bi foil, shown 
in Figure ^a). At the inflection point indicated by the arrow, /gj has its largest magnitude 
of —24.7 electrons. Uncertainties in /gj may arise from inaccuracies in the calculation, 
incorrect establishment of the incident energy, and the energy resolution. To account for 
these possibilities, the analysis was performed for deviations in /gj of about 20% (i.e., /gj = 
— 19.7 and —29.7 electrons). The results were incorporated into the error ranges tabulated 
below. 

V. X-RAY REFLECTIVITY DATA 

Reflectivity data at both energies are shown in Figure ^b) (symbols). These results were 
found to be both reproducible in energy and stable over time by measuring two 12 keV 
data sets, prior to and following the 13.4 keV measurements. The two 12 keV data sets are 
shown together as open circles in Figure ^b). The interference peak at = 2.1 is due 
to stratification of the atoms in planes parallel to the surface, a well established feature of 
liquid metals.— For < 1.8 A~^, the data exhibit a modulation indicative of a structural 
periodicity roughly twice that of the surface layering. This feature is consistent with Bi-In 
dimers oriented along the surface normal, which could also give rise to alternating Bi and In 
layers (bilayers) at the surface. We also find that the low-g^ reflectivity is strongly decreased 
when measured at the Bi Lm edge. The reflectivity decrease itself varies smoothly with 

and does not exhibit a bilayer-type modulation, instead suggesting that the surface Bi 
concentration is larger than that of the bulk. 

VI. ELECTRON DENSITY MODEL 

To investigate these possibilities, we calculate the reflectivity of a model density profile 
according to Eq. Q which is refined simultaneously against the data taken at both energies. 
Resonant effects are included by combining the model structure and the scattering amplitude 
into an effective electron density profile pcs{z). The energy dependence enters the analysis 
through the Bi concentration defined in pefr(-2), and also through the Fresnel reflectivity Rp, 
which is a function of the energy dependent mass absorption coefficient fi^^ and the effective 
bulk electron density poo oc {Z — /') that defines the critical angle qc- Table HI shows the 
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Bi22ln78: 12.5 74.5 7.66 17.5 — — 1.874 0.05154 

13.4 92.6 14.1 — — 1.764 0.05000 

"C. H. MacGillavry and G. D. Rieck, eds., International Tables for X-ray Crystallography Vol. 

Ill, Kynoch Press, Birmingham, England (1962). ^R. C. Weast, ed., CRC Handbook of Chemistry 

and Physics (see Ref. 20). "^B. L. Henke, E. M. Gullikson, and J. C. Davis, Atomic Data and 

Nuclear Tables 54 (1993) 181. 

TABLE I: Parameters used to calculate energy dependent x-ray reflectivity. 

values of these quantities used in our models. Reflectivity data acquired at each energy were 
normalized to the appropriate energy dependent Fresnel function. 

Following past practice,— our model incorporates layers of atoms having a Gaussian 
distribution of displacements from idealized positions nd along the surface normal direction: 

OO ^ 

Peaiz) = pooy^Fn 7=exp [-{z - ndf / al] . (2) 

The roughness arises from both static and dynamic contributions: 

2 2 2 _ / Q'max 



ires 



a^ = na + + ^ In ) • (3) 



Here a and (Tq are related to the surface layering coherence length and the amplitude of 
density oscillations at the surface. The last term accounts for height fluctuations produced 
by capillary waves, and depends on the temperature T = 101°C, the surface tension 7 = 
0.50 N/m,^^9 and wavevector cutoffs gmax = 0.99 A"-*^ and g^es ~ 0.024 A^-*^ (a slowly varying 
function of g^), as detailed elsewhere.— 

The scattering amplitude of each layer depends on the form factor and the effective 
electron density in each layer relative to the bulk, dependent on energy and Bi concentration: 

p Xn,Bi [fBijqz) + frnjE)] + (1 - Xn,Bi) [flnjqz) + fU a^n.BiPoo.Bi + (1 - a:n,Bi)Poo,In 

Poo, bulk 

(4) 
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FIG. 2: (a) Model surface- normal number density profiles, relative to bulk atomic percent. ( — ): 

total number density; ( ): In density; ( ): Bi density, (b) Calculated Presnel-normalized 

reflectivity curves. ( — ): 12.5 keV; ( ): 13.4 keV. i: Surface layering with uniform composition 

(xBi = 0.22). ii: Surface layering, with Bi enhancement in the first atomic layer (xo,Bi = 0.35). 

For the first few layers (n=0, 1, 2), the weight Wn may differ from unity and the Bi fraction 
3^n,Bi can vary from the bulk value = 0.22. 

VII. DATA FITTING ANALYSIS 



We now describe the ingredients that are required to fit the data with this model The 
interference peak at = 2.1 A^^ can be reproduced by a simple layered profile in which 
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Figure d W ao wq xq^bi wi xi^Bi ^2 X2,Bi 

1 2.81 0.54 0.64 0.98 0.35 1.01 0.22 0.98 0.23 

2i 2.81 0.48 0.64 1.0 0.22 1.0 0.22 1.0 0.22 

2ii 2.81 0.54 0.64 1.0 0.35 1.0 0.22 1.0 0.22 

TABLE II: Fit parameters for model profiles, identified by the figure in which they appear. The 
length scales d, a, and ctq are in units of A. 

Xn,Bi = XBi for all n. The Bi and In number densities for such a model are shown in 
Figure |2fa)i, along with their sum, the total atomic fraction relative to the bulk. The cor- 
responding reflectivity curves calculated for both x-ray energies are compared to the experi- 
mental data in Figure Efb)!. Since the Bi concentration is uniform, the energy dependence is 
so slight that the curves overlap almost completely on the scale of the figure. Turning again 
to the data, the reduced intensity in the region < 1.8 A^^ when measured at the Bi Ljjj 
edge indicates a reduction in the scattering amplitude at the surface. This implies that the 
Bi concentration is enhanced there. Increasing the Bi fraction from 22 at% to ~ 35 at% in 
the first layer (n = 0) produces an appropriate energy dependence (Figs. |2ta)ii and|21^b)ii). 

Although at this point the frequency of the low-g^ modulation is not well described, 
the fit is considerably improved by allowing the Bi fraction and total number densities to 
vary for the first three surface layers (n=0,l,2). We find that the Bi fraction for n=l,2 is 
essentially equal to the bulk value of 22 at%, while the total number densities for n=0,l,2 
have values of 0.98, 1.01, and 0.98, respectively (Figs. [Tl^b),(c)). Thus, the detailed shape 
of the low-g^ modulation is modelled by a very slight density wave of about 2% affecting 
the amplitudes of the first few surface layers. Increasing the number of model parameters 
to allow for shifts in positions and widths of the surface layers resulted in marginally better 
fits, but at the expense of high frequency Fourier components appearing as small wiggles 
in the calculated reflectivity. Variations in these extra fit parameters are extremely slight, 
and we doubt whether they have any physical basis. Parameters for all models are shown 
in Table HH 

This analysis demonstrates that to model the essential features of our data, there is no 
need to invoke long-range compositional ordering on a second length scale in the surface- 
normal direction, which we had suggested based on the previous non-resonant reflectivity 
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measurements.— Still, we thought it important to investigate additional, specific models 
based on Bi-In pairs oriented along the surface normal. To test for alternating Bi-rich 
and In-rich layers, we attempted fits in which we forced the Bi and In compositions to be 
substantially different from the results shown above. We also described Bi-In pairing by 
allowing the positions, but not the densities, of the surface layers to vary. None of these 
profiles successfully described the data. 



VIII. SUMMARY 



Our principal finding is the Bi enrichment of 35 at% in the surface layer, compared 
to the bulk value of 22 at%. This is considerably less Bi than would be expected in the 
absence of attractive Bi-In interactions, which can be estimated from the surface free energy 
requirement)^ 

7in + In - - XBi - = 7Bi + In " 7 ^ - a^Bi - • 5 

a \ 1 — / 4 a a \ xbi / 4 a 

The quantity w is the excess interaction energy of Bi-In pairs over the average of the Bi-Bi 
and In-In interaction energies; for an ideal mixture, w = 0. This analysis assumes that the 
inhomogeneous region is confined to a single atomic layer, the atoms are close-packed and 
take up an area a, and w is small. Extrapolating the measured surface tensions to 100°C,— 
7in = 0.56 N/m and 7Bi = 0.41 N/m. Using the Bi atomic size, a = 7r(3.34/2)2 (for In, 
the atomic diameter is 3.14 A),— and the bulk composition xbi = 0.22. The equilibrium 
surface composition xq.bi is found from the intersection of plots of both sides of Eq. El For 
w = 0, this analysis predicts a surface segregation of 69 at% Bi (Figure |21^d), dashed lines). 
To reproduce our experimental finding that xo,Bi = 35 at%, w must be negative, with a 
magnitude of ~ lOfc^T (Figure Efd), solid lines). Although this large value of w is most 
likely outside the range of validity of Eq. the analysis certainly illustrates the qualitative 
effect of attractive heteroatomic interactions on the surface composition. In this Bi-In alloy, 
pairing does in fact defeat Gibbs adsorption in the sense that the surface energy is optimized 
not by segregating a large fraction of Bi, but by forming larger numbers of Bi-In neighbors in 
the surface layer. Exactly how this balance plays out in Bi-In alloys with the stoichiometric 
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bulk compositions Biln and Biln2 remains to be seen. 
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